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. . . We derive an effective polaronic interaction Hamiltonian, exact to second order in perturbation, 

' for the spinless one-dimensional Holstein model. The small parameter is given by the ratio of the 

, hopping term (t) to the polaronic energy (g^LOo) in all the region of validity for our perturbation; 

however, the exception being the regime of extreme anti-adiabaticity {t/u>o < 0.1) and small electron- 
phonon coupling {g < 1) where the small parameter is t/uiQ. We map our polaronic Hamiltonian 
onto a next-to-nearest-neighbor interaction anisotropic Heisenberg spin model. By studying the 
mass gap and the power-law exponent of the spin-spin correlation function for our Heisenberg spin 
model, we analyze the Luttinger liquid to charge-density-wave transition at half-filling in the effective 
polaronic Hamiltonian. We calculate the structure factor at all fillings and find that the spin-spin 
correlation length decreases as one deviates from half-filling. We also extend our derivation of 
polaronic Hamiltonian to d-dimensions. 
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I. INTRODUCTION 



Understanding the many-body aspects of the Holstein model^ has been a long standing open problem. Significant 
progress has been made in understanding the single polaron problem through analytic treatments in the small and 
large polaron limits and numerical approaches in the intermediate size case^. While studies of the many-polaron 
Q problem, involving spin degrees of freedom, yielded interesting insights for bipolarons and their phase transitions^'"^, 
O ' the simpler case of spinless many-polaron problem has received only scant attention^. With the renewed interest in 
strongly correlated manganite systems^ (which are of high electronic density when the electron-phonon interactions 
' are supposed to be important) it is imperative that the effective interaction between polarons be understood so that 
^ , a serious attempt at explaining the rich phase diagram of these systems be made. Although studying manganites 
• demands knowledge of the effective Jahn- Teller polaronic interaction, in the low-doped case the effective Hamiltonian 
T— I ■ at K for the occupied states can be taken as a Holstein model^. Thus a good understanding, involving exact results, 
'/"J of the simpler effective polaronic interaction for the Holstein model, which has been elusive so far, is highly desirable. 
Furthermore, an effective polaronic Hamiltonian even in the simplest spinless ID case would also be quite useful for 
modelling Luttinger liquid (LL) to charge-density-wave (CDW) transitions in half-filled systems. Quasi-ID organic 
conjugated polymers [such as (CH)x], charge transfer salts [such as TTF(TCNQ)], and inorganic blue bronzes (e.g., 
Ko.sMoOs)^ are good candidates for such broken symmetry in the ground state leading to unit cell doubling. 
' The present paper, using Lang-Firsov transformation* , provides a transparent perturbative approach to deriving the 
^ ' effective Hamiltonian of interacting polarons in d- dimensions when the band narrowing is significant. The resulting 
I ■ polaronic Hamiltonian, exact to second order in perturbation, is studied in ID at K for density-density correlation 
' O ] effects. The correlation function exponent of the concomitant quasi-long range order is demonstrated to be useful in 
^ I characterizing the Luttinger liquid and the charge-density-wave phases of the system. 
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II. EFFECTIVE POLARONIC HAMILTONIAN 



^ ' We begin in ID by taking the unperturbed Hamiltonian to be the non-interacting polaronic term^ 
^ ! 

j 3 

-Ji^(c]c,+i+H.c.), (1) 

3 

with the perturbation being 

H' = Y. H3 = -Ji Y.^c]c,+,{S'^Sl - 1} + H.c), (2) 
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where Hq + H' make up the Lang-Firsov transformed Holstein Hamiltonian and Cj (aj) is the fcrmionic (phononic) 
destruction operator, wq the Debye frequency, Sj_ = exp[±5(aj — Oj+i)], Ji = texp{—g^), with t being the hopping 
term, and g'^iVo the polaronic binding energy. The eigen states are given by |n,m) = \n)ei (8 |m)p/i with |0,0) being 
the ground state with zero phonons. The cigcn energies arc £"„.,„ = E^j + EfJ^ . Since (0, 0) = 0, the first 

order perturbation term is zero and the relevant excited states correspond to states with non-zero phonons. Next, we 
represent \m)ph in real space as phononic excitations at different sites with one phonon state being a]|0)pf, which is 
A^-fold degenerate and can correspond to any site j. On the other hand, the electronic state is represented in 
the momentum space. We will now calculate second-order perturbation terms exactly 

Now E*^ —Eq ~ Ji and ^E^ = E^ — E"^^ is a non-zero integral multiple of ujq. Assuming Ji « ujq (which certainly 
is true for realistic values of 2 < t/wo < 4 and 6 < < 10 found in manganites) and using \'n){n\ = I we get the 
corresponding second-order term in the effective Hamiltonian for polarons to be 

= ^ ^ {0\phHi\m}ph{m\phHj\0}ph 

Lj m 

In the above equation, the term Hj produces phonons at sites j and j -|- 1. Hence to match that its counterpart Hi 
should produce phonons in at least one of the two sites j and j + 1. Thus the index / = j — 1, j, or j + 1. Next on 
defining P±{j;m) = {0\phS± - 1 1^)^/1 and bj = cjcj+i we get 

H^'^ =-y2^[{b%P+U;m)+b,.,b,P.U - l;m) 



+bj+ibjP-ij + 1; m)}P[(i; to) 
+{bjb]P^{j; to) + b]_,b]P+{j - 1; to) 

+6]+i6]P+(j + 1; TO)}P!(j; to)]. (3) 

Then using (a'f)"|0) = \/n!|n') with \n') being a state with n phonons we get the effective polaronic Hamiltonian to 

be 

^e// = -g^'^o J2 ^1 ~ "^1 + H.c.) 

j 3 

+r ^ UjUj+x + 2 Ja ^{c]_-^njCj+x + H.c.) 

3 3 

-J2 ^(4-1^^+1 + H-c.) - E"^' (4) 

J i 

where ^ ^^m9) + V2{a)l and ^ ^J^{g) with h{g) ^ E^^i £ ^nd h{g) ^ EZi E™=i J^nZ) - 
is of interest to note that the single polaronic energy part of the above Hamiltonian matches with the self-energy 
expression at = obtained by Marsiglio^" and the self-energy at a general k by Stcphan^^ and lends credibility to 
our results. Furthermore, in the work of Hirsch and Fradkin^^, the coefficient of nearest- neighbor interaction agrees 
with our coefficient for large values of g while the coefficients of the next-to-nearest-neighbor (NNN) hopping 
are in disagreement with those of ours and the results of Refs.^°'^^. Next we make the connection that, on using 
Wigner-Jordan transformation = nj<j(l — 2nj)ct, we can map the effective polaronic Hamiltonian exactly (up to 
a constant) on to the following NNN anisotropic Heisenberg spin chain: 

= -a'^o E -I - ^1 E(->j;i + H-c.) 

3 3 

E ^M+i - E«i^7+i + (5) 
j j 

where the coefficient of the first term represents coupling to a longitudinal magnetic field. Although the NNN 

interactions in the above Hamiltonian do not produce frustration, nevertheless the Hamiltonian cannot be solved by 
coordinate Bethe ansatz^^. Hence we take recourse to analyzing the properties of the effective Hamiltonian numerically 
by using modified Lanczos technique (see Ref.^^ for details). 
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III. LUTTINGER LIQUID TO CDW TRANSITION 



The spin Hamiltonian H^J^P in Eq. (5) with J2 = 0, i.e., without NNN interaction, has been shown to undergo a 
LL to CDW state transition at zero magnetization (t| = 0) when = 2 Ji and at non-zero magnetizations is 

always a LL^^. On including a non-zero J2, the disordering effect increases because the NNN interaction is only in 
the transverse direction and the LL to CDW transition will occur at higher values of J^. We expect that including 
J2 does not change the universality class and that the central charge c = 1. 

We first study the static spin-spin correlation function on rings with even number of sites N and extract information 
about the critical exponent. The static spin-spin correlation function for a chain of length N is given by Wi{N) = 

TfY^f=\{SiSi+i) ^nd has the asymptotic behavior X\m.j^^ooWi{N) rs ^''pP for the anisotropic Heisenberg model 

when I ^ 1"'^^. Furthermore, A is an unknown constant and 1 < 77 < 2 when the system is in disordered (or LL) state, 
= 1 is the transition point to antiferromagnetic (CDW) state, and 77 = means system is totally antiferromagnetic 
(or CDW). 

We will now derive an analytic expression for the critical exponent rj based on the work of Luther and Peschel^^. 
The effective polaronic Hamiltonian given by Eq. (4) can be written in momentum space as 



TTpol 

^eff 



-2 Ji J2 cos{k)clck + ^ 51 t^os(g)p(g)p(-g) 



4J2 
N 



^ C0s(fc + fc')4+g4,Cft/+gCfc 



-2J2 ^ cos(2fc)4cfe, 



(6) 



X]p c^^pCp. Furthermore, constant terms have been ignored. Next, we linearize the kinetic energy term 



where p{q) = 4 

close to the Fermi points and follow it up with the bosonization procedure. Then on taking exchange effects into 
account, as pointed out by Fowler^^, we obtain the following bosonized Hamiltonian 



TTpol 

^bos 



477 Ji + 4 + 8 J2 



+ 



N 

8 - 32 J2 
N 



XI Pi{k)pi{-k) 



fe>0;i=l,2 



Y,Pi{k)P2{-k). 



(7) 



fc>0 



It is important to note that only the forward scattering part involving the coefficient contributes to the self energy 
correction. 

Now, to calculate the critical exponent r?, we will follow the usual procedure and diagonalize the bosonized Hamil- 
tonian of Eq. (7) using the following transformations. 



Px (k) = pi{k) cosh + p2 {k) sinh (p, 



and 



P2{k) = P2{k) cosh^ -|- pi(fc) svahcj). 
Then, on setting the coefficient of the off-diagonal term equal to zero in the transformed Hamiltonian, we get 

2 - 8J2 



tanh 2(f) 



2nJi + 2 + 4 J2 ' 



(8) 



Using Eq.(8) we obtain 



2e^'^ = 2 



1 + 



1+ 2J1_ 2i 



(9) 
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IV. RESULTS AND DISCUSSION 



It is known that for an anisotropic Heisenberg spin chain, when N/2 is even the correlation function goes to zero 

smoothly as the longitudinal interaction goes to zcro-'^*. Hence we have calculated T4^Ar/2(iV) only for odd values of 
N/2 with = 6, 10, 14, 18, and 22 at Ji = 1 and different values of and J2- Using a linear least squares fit for a 
plot of log W;v/2(-^) versus log A'' we obtained the value of r] from the slope at each value of and J2. The error in 
77 for all curves is within « 0.05 and hence verifies that one has the expected quasi-long range order. For = 2 and 
J2 = we get rj = 0.96 ± 0.05 which includes the exact value of 77 = 1 obtained from Bethe ansatz. Thus we expect 
the T] values obtained by our procedure to be reasonably accurate. Since for > 2 and J2 = we obtain a CDW 
state, by increasing J2 at any > 2 we should increase the disordering effect and hence we see in Fig. 1 that r] value 
increases. We find that for « 6 the value of 77 becomes slightly negative but with magnitude within the error of 
0.05. At higher values (w 10 and higher) r/ tends to zero. At small values of (< 0.5) as J2 increases initially 77 
increases even to values above 2 and then decreases to values below 2. We think that this interesting feature is due to 
smaller values of J2 enhancing the disordering effect while the larger values of J2 increase correlations build up with 
the system becoming less LL like. However the behavior at = and J2 > needs further imderstanding. Our 
derived analytic expression, reliable at small values of Ji/ J\ and rj < 2, shows that 77 does increase with increasing 
values of J2 and gives values reasonably close to the numerical ones for J^/Ji < 1. 

We will now consider the mass gap, at the half-filled state for the Hamiltonian H = H^j^^ + g'^ujQ Uj [see Eq. 
(4)] , defined as twice the energy difference between the ground state with 1 + N/2 electrons and the ground state of 
the N/2 electronic system. The mass gap is calculated for rings with N = 10, 12, 14, 16, 18, and 20 sites. Including N 
= 6 and 8 only increases the error. The mass gap plot in Fig. 2 is obtained using finite-size scaling by plotting mass 
gap versus 1/N and extrapolating the linear least square fit to the value corresponding to 1/N = 0. In the plot the 
size of the symbol is larger than the error. From the inset of the plot of mass gap versus at various values of J2, we 
see that the mass gap goes to zero at « 1.4 at J2 = which is a significant underestimation of the transition value 
of = 2. Also on comparing with Fig. 1, we again notice that the LL to CDW transition value of at different 
J2 values is grossly underestimated. Furthermore, as expected, mass gap increases (decreases) monotonically with 
(J2) at a fixed J2 (J^). 

We will now discuss the region of validity for our model, given by Eq. (4) and as depicted in Fig. 3 (region above 
the lower curve), in the two-dimensional parameter spac;e of g and t/uJo- Firstly, since we use the assumption that 
Wo ^ Ji in our derivation, we choose the validity condition as loq > 10 Ji. Next, we would like the second order 
energy term E^"^^ in the perturbation series to be much smaller than the unperturbed term Eo^. We find that for 
t/^o ^ li the condition = 10 Ji produces a boundary on which the ratio Eq^/E^"^^ > 5 with the ratio increasing 
rapidly as t/uJo decreases. As for t/uJo > 2, we find that the condition g^LOo > 3J^ is more restrictive than the first 
one {u>o > 10 Ji) and produces a ratio of Eq^q/E^^^ > 3 at the boundary. Next we will discuss the LL to CDW phase 
transition boundary obtained from 77 = 1 condition and depicted by the upper curve in Fig. 3. We find that the 
phase transition points lie within the region of validity only for t/coo < 0.6. In the region to the right of the dashed 
vertical line and below the region-of-validity curve the phase boundary cannot be determined using our model. It is 
important to note that the experimentally realistic pa,ram,eter regime 6 < < 10 a,nd 2 < t/uiQ < 4 lies mostly inside 
the region of validity. Upon comparing our numerical phase transition results with those of Bursill et al}^ , we find 
that for small values oit/coo < 0.1 the critical gc values agree well. However for larger values the results do not agree. 
At t/uo = 0.5 our gc = 1.45 ± 0.02 (with Eq^o/E^"^^ > 17 and w/Ji > 15) is noticeably smaller than the gc = 1.63(1) 
of Ref.-^® and at t/ujo = 1 we find that gc < 1-52 whereas Bursill et al. get gc = 1.61(1). As for t/cjo > 2, our region 
of validity lies above the phase transition boundary given in Ref.^''. However, interestingly, the numerical estimates 
of the critical gc by Hirsch and Fradkin^^ are consistent with our results with their gc value agreeing with ours at 
t/uJo = 0.5, while at higher values oit/uiQ their gc values lie outside our region of validity. 

Now that the region of validity has been identified, we will analyze within this region the small parameter for 
our pertubation theory. For g > 1, one approximates fi{g) ~ eP / and [2/1(5) + !i{9)\ ~ I '^9^ with the 
approximations becoming exact in the limit g 00. Then the effective polaronic Hamiltonian of Eq. (4), for the case 
g > 1, simplifies to 

H^^ff ~ -9'u;o[J2 + Ce"^' Y.^c]c^+^ + H.c.) 

3 3 

+C'e-''' E{4-i(l - 2n,)c,+i + H.c.} 

3 

3 
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where C = t/g^ujQ is the polaron size parameter. In the region of vahdity for our model, when the adiabaticity 
parameter t/wp > 0.2, we have the constraints g > 1 and g^uji^ > 2t. Thus we see that, for the region t/ojQ > 0.2, the 
polaron size parameter C is the small parameter. 

Now, for the extreme anti-adiabatic regime of t/u)o < 0.1, all values of g are allowed by our model. When g > 1, 
again Eq. (10) is valid with the same small parameter (. However, when g < 1, we make the approximations 
hid) ~ 9^ and [2/1(5) + /2(g)] ~ 2/1 (gi) with the approximations becoming exact in the limit 5 — > 0. Then, for g <1, 
the effective polaronic Hamiltonian given by Eq. (4) becomes 

H^^/f ~ -g'uJoiY^ + Ce-^' ECS'^^+i + H-c) 

i j 

+ (^) ' ''''' E{S--i(l - 2n,)c,-+i + H.C.} 

+4(-iye-2s=^n,(l-n,+i)]- (11) 

Thus, for the regime t/uo < 0.1 and g < 1, the adiabaticity parameter t/wQ is the small parameter in Eq. (11) with 
5^1 corresponding to small polarons and 5 — > (such that g'^too «t) corresponding to large polarons. 

Finally, we also study the static structure factor Sis[{k) = '^^^i e^''^\Vi{N) . The structure factor offers information 
about the correlation lengths even in LL phase at all filling factors n. In fact, the correlation length decreases 
with increasing width of the structure factor near its peak at 27rn. We first observe that X)fe'^Af(^) ~ ^ ^^'^ ^^^^ 
Sn{0) — 47V(n — 0.5)^ all of which arc borne out by both the plots in Fig. 4 done at = 20. The plots arc only for 
k = 2mn/N with m = 0, 1, ...,N/2 as Swik) is symmetric about n and are only for n < 0.5 because of particle-hole 
symmetry. Fig. 4(b), plotted for t/uo = 0.5 and g = gc = 1-45 (or = 2.53 and J2 = 0.245), corresponds to LL- 
CDW transition point at n = 0.5, while Fig. 4(a), done for the realistic values of t/ujo = 3 and g = 3 (or « 3000 
and J2 ~ 0.38), depicts situation deep inside the CDW phase at n = 0.5. Now, we know that at n = 0.5, the structure 
factor Sn{tt) ~ / and hence diverges for 77 < 1 (CDW regime) as iV 00 with the divergence being faster as we go 
deeper inside the CDW regime. Also the structure factor remains finite at 27rn for all other filling factors even when 
N ^ 00 because here the system is always a LL. From plot (a) we infer that deep inside the CDW state S'Af(7r) « TV 
and SN{k 7^ tt) « 0. As for the CDW transition point depicted in plot (b) at n = 0.5, the structure factor peaks 
sharply but more gradually at fc = tt. In both plots the largest peak occurs for n = 0.5 with peak size diminishing 
and curve width increasing as values of n decrease. Thus, we see that for n ^ 0.5 also, short range correlations exist 
with correlation length decreasing as deviation from half-filling increases. 

Lastly and importantly, using arguments similar to the ID case, we have also derived the effective polaronic 
Hamiltonian in d-dimensions to be^*^ 

-0.5J^^n,(l-n,+,), (12) 
j,s 

where d corresponds to nearest neighbor^"'^. 

In conclusion, we have derived an effective polaronic Hamiltonian for the spinless ID Holstein model which is found 
to be valid in most of the experimentally realistic regime. We mapped the effective Hamiltonian onto a next-to- 
nearest-neighbor anisotropic Heisenberg Hamiltonian. Using modified Lanczos technique extensively, we computed 
the static spin-spin correlation exponent 77 and the mass gap at half-filling for general values of the parameters in 
the effective spin Hamiltonian. The mass gap values were found to significantly underestimate the critical electron- 
phonon coupling gc- In contrast, the r] values were found to give reliable estimates of gc and consequently were used to 
determine the LL-CDW quantum phase transition. The structure factor calculations revealed that correlation length 
diminishes with increasing deviation from half filling. Lastly, our approach, exact to second order in perturbation, 
is extended to obtain an effective polaronic Hamiltonian in d-dimensions also. Our work opens up a whole host of 
future challenges such as: (1) Extension to finite temperatures and studying metal-insulator transition; (2) Including 
Hubbard interaction U; (3) Analyzing the d-dimensional model in Eq. (12)^''; and (4) Deriving analogous effective 
Hamiltonian for Jahn- Teller systems. 
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FIG. 1. Plot of spin-spin correlation exponent for various values of and J2. The dashed lines are guides to the eye. 
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FIG. 3. Plot of region of validity boundary and LL-CDW phase boundary for g versus t/uio- The errors are smaller than the 
symbols. The crosses depict realistic regime. 
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FIG. 4. Structure factor plots at various values of k and n. 
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